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ABSTRACT: Whether for industrial, transportation or complex systems, ball bearing are considered as major 
factor in kinematics rotating systems. Apart from their strategic role in the drive assemblies, these mechanical 
components can require an increased rate of reliability. Bearings turbines and rotating machines are primordial 
in most mechanical components which contribute directly to the performance of automotive and aerospace 
engines through reliability. Due to their role as liaison between fixed and moving parts, any failure could have 
catastrophic consequences such as loss of engine use. This research is aimed to develop the analytical 
fundamentals details of new model which is dedicated to the dynamic behavior of rotating ball bearings. 
Theoretical research will be detailed enough to provide a simple model that lends easily to programming, but 
powerful enough to incorporate the effect of a maximum of influential parameters on the vibration of 
bearings. The essence of this work is to present an analytical model of a specific type single of row effect balls 
(SKF 6004). This study is a first step dedicated to define the geometric characteristics, and in second time to 
determine the equations governing the distribution of responsibilities within the ball bearing. 
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I. INTRODUCTION 

Even if the geometry of a ball bearing is perfect, it will still produce vibrations. The vibrations are 
caused by the rotation of a finite number of loaded rolling contacts between the balls and the guiding rings. 
Because these contacts are elastic, the bearing stiffness becomes explicitly dependent on time. In general, a time 
varying stiffness causes vibrations, even in the absence of external loads. Since the stiffness can be regarded as a 
system parameter, the variable stiffness leads to a so-called parametric excitation. It is one of the major sources 
of vibration in ball bearings. The first systematic research on this subject was conducted by (Perret, 1950) and 
(Meldau, 1951). 

(Dowson and Higginson, 1977) have proposed the first analytical solutions to EHD in case of 
cylinder/plane contact. Then, in the middle of 1970s, digital solutions contact ellipsoid/ maps have been 
proposed by (Hamrock and Dowson, 1977). Thus, using these procedures, we can go back to the thicknesses of 
lubricant in contact and we can realize that these thicknesses vary between a few nanometers. From the 
viewpoint of contact fatigue in the presence of indentation, (Wang and al., 2007) exhibit a work based on the 
introduction of a ceramic rolling element (Si3N 4 ) from the rolling ones in steel. This study demonstrates a 
"smoothing effect" created by the passage of ceramic rolling element which generates strong plastic 
deformation of bead. This "smoothing effect" makes it possible to reduce the height of beads and, at the same 
time, local overpressure. This suggests that the addition of ceramic allows increasing number of cycles to 
chipping adjacent indents. In addition, the work presented by (Jacq, 2001) shows a smaller effect of the load on 
damage by chipping in vicinity of the with respect to sliding. (Cheng, 1997) present a model including contact 
fatigue crack initiation in slip bands formed in the grains of the material. The slip bands are modeled by two 
stacks of dislocations of opposite signs, of which the length is equal to that of the grains. Accumulation of 
dislocations occurs during the cyclic loads due to irreversibility of dislocations movement. This type of 
analytical study is actually based on the work of (Tanaka and Mura, 1981), together with its collaborators 
studied crack initiation in a crystalline material subjected to cyclic loading. Moreover, the modeling of the 
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propagation of cracks in contact fatigue was studied by many other authors. Of particular we note: (Keer and al., 
1983), (Lee and al., 1986), (Murakami and al., 1994). 

II. ROLLER BEARING 
The bearing is a mechanical member permitting rotation of one part relative to another by reducing 
friction. It has three main functions: 

• -The transmission of forces 

• -Rotation states in one piece with minimum friction 

• -The position of a shaft relative to the housing 

• Thus, if we describe a ball, we find rolling elements, separated by a cage, an inner ring and an outer ring 
provided with raceways (Figure 1). 

The ball bearing can replace the sliding friction by rolling one. However, various body contacting (inner ring 
and balls for example) showing differences speed in their contact areas that will generate local slipping. Indeed, 
the friction generated by the fastest surface will tend to pull the slow surface. 




outer ring 



Figure 1 : Overview of ball bearings components 




To ensure a long life, the roll bearings are lubricated with grease or oil. According to their operating 
conditions (load, speed ), different regimes of lubrication can be occurs: 

- A system of hydrodynamic lubrication (HD) 

- A system of elastohydrodynamic lubrication (EHD) 



2.1. The E.H.D Contacts In Bearing 

A roll bearing can force the position of shaft relative to system. The rolling elements transmit the 
forces applied by the shaft to the inner race fixed to the outer ring. They facilitate also the rotation of shaft 
relative to housing (Figure 2). 




Figure 2: Schematic operation of a roll bearing 



Geometry of rolling elements and raceways (inner and outer ring) differs depending on the loading conditions 
(axial and / or radial moments) and speed. Bearings angular contacts are designed to simultaneously collect 
radial and axial loads (Figure 3). They can operate at very high speed and are therefore used to stopping gas 
turbines. The axial component of the thrust and the radial component generated by gas is due to the own weight 
of the shaft and the unbalance effects engagement in the case of the use of a speed reducer. 
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Figure 3: Ball bearing angular contact 



Because of the high speeds (above 20000 trs/min ) , centrifugal forces (Fc ) acting on the balls have the 
same order of magnitude, as the forces associated with contact with the rings ( Figure 4). These efforts are 
reflected in resolution of load distribution in the bearing, as the gyroscopic moments acting on the ball. The 
details of the corresponding quasi- static analysis are described by (Harris, 1984).The determination of bearing 
kinematics is based on assumption of pure rolling, which is defined by the equality of geometric points of speed 
ball and ring contact considering non-deformable solids. The spin rate of the balls (co b ) and the orbital velocity 
(co m ) depend on angles of contact between ball and the outer ring (a D ) and ball with inner ring (a;) (Figure 4). 
Centrifugal forces tend to flatten ball to the bottom of outer ring and corresponding contact angle (a 0 ) is less 
than the angle of geometric contact (a 0 ). The opposite occurs for the inner contact angle (a ; ). The pivot bearing 
is impossible without simultaneously ball contacts / inner ring and ball / outer ring. Due to centrifugal forces, 
pressures and contact surfaces are more important to touch ball / ball at the outer ring / inner ring contact. The 
pivot is then assumed to be zero, thus fixing the kinematics of the bearing. This is assuming control of the outer 
ring ball (Dawson and al. 1986). 




Figure 4: Distribution of kinematic and load bearing ball angular contact with high speed 

The contact between the rolling elements and rings is characterized by non-compliant type ball / plane, 
where the local pressure has the effect of deforming the solid surfaces. Hertz has been the pioneer in this field 
(1881), laying the foundations of the Hertzian contact theory which estimates the distribution of pressure, 
deformation and dimension of contact. The speed differences between the solids in contact, due to sliding of 
pivot, are responsible of frictional forces, directly proportional to the coefficient of friction. This reflects the 
relationship between normal and transmitted tangential at contact forces. In the case of a dry contact, the friction 
coefficient is higher than coefficient of lubricated contact. It leads to friction levels which is too important in roll 
bearing and cause more or less long term ruin in operation. 

2.2. The H.D. Contacts In The Bearing 

In a ball bearing exposed to axial and radial forces, the charge distribution is not in cyclic symmetry. 
The contact angles are different for each of rolling elements which have not all the same kinematics. However, it 
is necessary that all rolling elements retain on average a regular spacing from each other to ensure the proper 
functioning of mechanism service. The cage can meet this need by positioning each element in a cell. It supports 
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the efforts made by opponents of the printed contacts with rings (Figure 5) movement. The lubricant can 
accommodate the difference in speed between the rolling element and the cage and avoid a dry friction between 
the contacting surfaces. 



1 

y Cage pocket 




1 t 

Figure 5: Description of contact ball cage in a ball bearing 

The cage plays the role of separator and must, however, be guided in rotation especially in high speed 
applications where the imbalance caused by the eccentricity of gravity center can be fatal. It is possible to guide 
the cage by the inner ring or the outer ring (Figure 6). In the first case, the cage will tend to be driven (mobile 
ring) but is sensitive to unbalance and therefore requires a good balance. In the second case the cage is 
restrained by the outer ring (fixed) and is more sensitive to the phenomenon of torque (lash adjuster operation 
by differential centrifugation and expansion). 



Inner ring 

1 Outer ring 




Figure 6 - Guiding of the cage and inner ring/outer ring 

III. MATHEMATICAL DEVELOPMENT OF A SINGLE BEARING 
3.1. INITIAL MODEL 

In this part we consider a rotor, bearing rings and support. The model becomes complex with a large 
number of degrees of freedom (ddl). In our model, we applied three degree of freedom system as described in 
Figure 8, in order to focus the study on the vibration response of the balls (Sanchette-Gosset, 1993) and (Tanaka 
and Mura, 1981). The dynamic behavior of rigid bearing support and the contribution of the support is also 
considered. Thus, this simplified model that will be adjusted with the calculation of vibration response. 



68 



Analytical Approach And Mechanical Estimate 




Outer ring Fluid film BaU Fluid film Inner ring 



M 




Figure 8: Model of bearing system 

With: 

Ml : Mass of outer ring [kg] 
M3: Mass of inner ring [Kg] 
M2: Mass of ball [kg] 
Kl: stiffness of outer ring [N/m], 
K3: stiffness of inner ring [N / m] 
K21: Stiff ball [N/m] 

C22: Damping ratio Film external fluid [N.s / m] 

C21: Damping coefficient of internal fluid film [N.s / m] 

The different parameters governing the behavior of this model were developed using the following analytical 
elements: 



The Rings :We considered a mass-spring system with neglected structural damping, the rigidities of the two 
rings is determined by the following equations: 



Rigidity of the inner ring 



(1) 



Rigidity of the outer ring 1 1 1 (2) 

Each of inner and outer rings will be represented by a mass-spring combination which corresponds to a single 
degree of freedom system. The internal damping has been neglected. The natural frequency of the flexural 
vibration mode number "n" is given by the following expression: 

n[n 2 



(o.. = 



+ n 2 



•l] \m_ 

l juR 4 



(3) 



where: 

con angular frequency of mode [rad/s] 
n order mode 

E modulus longitudinal elasticity [N/m2] 

I moment of inertia cross section of the ring [m4] 

li linear density [kg/m] 
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R radius of ring [m] . 



Modes 0 and 1 are considered as rigid modes. Therefore the second mode is considered as the first bending 
vibration mode. The angular frequency of the system to a single degree of freedom will be taken as: 



C0„ = 2.68. 



— 

I MR* 



THE BALLS : Rigidity ball has the order of 8.3 109 N / m. Compared to other rigidities, it is considered 
infinitely rigid. Therefore modeling the ball was deemed to be a mass element. 

THE FLUID FILM : The development of elastohydrodynamic theory of lubrication (EHD) showed that films 
with a few micrometers thick are occurs in rolling contact. Therefore, the characterization of lubrication system 
depends on the value of lubricant film thickness. The theory developed by Hamrock [1], expressed in terms of 
stiffness and damping coefficients dimension film of fluid is given by the following expressions: 
Fluid film stiffness: 



k = - 



sin 2 ^, +- 



3xs n 



Amortization of fluid film : 



-sin d n cos L +- 



(l-^o) 3 



0/ cos^ 0 ] 



C = 



-[-, 
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2x3/2 



w r r k 2(i-^) 

With: 



sin 2 ^ + 



4£ n 



sin4cos4+— — —cos 4 



2d- sir 



s Q [16s 2 Q +x 2 (l-s 2 0 )] 



1/2 



cos 2 (fi 0 



2x1/2 



h 

c 



tan ^ 0 



s 0 =1 



And: 

Xk defines the ratio of width length of contact area 

(p0 is the maximal angle, described by the distribution charge 

Wr is the resultant force applied to fluid film 

£ 0 is an eccentricity coefficient 

c is the diametrical bearing clearance 
h is the thickness of the fluid film 



(5) 



(6) 



IV. RESULTS AND DISCUSSIONS 

The general equation to solve a system of three degrees of freedom subjected to three forces, depending 
on the mass matrices, damping and stiffness being have a general form: 

[A\.{y} + [Bl{y} + [C].{y} = {F} 

A, B and C represent respectively the masses of balls, damping and stiffness system. 



A = 



0 
0 



0 

M 2 
0 



0 
0 



B = 



r 

^22 


-C 

^22 


0 


-C 

^22 


c +c 

^22 T ^21 


-C 

^21 


0 


-C 

^21 


r 

^21 
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C = 



v 12 



/v 22 
^22 ^21 



0 



0 

~^21 
fc 3 + K 21 



The motion vector {y} and the force vector {f}. The vector F represents the impact forces generated by the fault 
of passage balls. 



M = 



^3 



3 J 



The equation of motion, governing this system model is the following: 

M iy\ +C 2 2 yi ~ c 22y 2 +( k i +^22)^1 -k 22 y 2 = F i 

M 2 y 2 -C 22 y 2 +(C 22 + C 2l )y 2 - C 2l y 3 -k 22 y x + (k 22 + k 2l )y 2 -k 2l y 3 = F 2 

M 3 y 3 +C 2l y 3 -C 2l y 2 +(k 3 +k 2l )y 3 -k 2l y 2 = F 3 
(7) 

This can be written in matrix form: 



M x 


0 


0 " 








" r 

^22 


-C 

^22 


0 








fcj + £ 22 


-k 

^22 


0 




r y^ 






0 


M 2 


0 


« 


y\ 


> + 


-C 

^22 


c +c 

^22 T ^21 


^21 


< 




> + 


^22 


k 22 + k 2i 


~^21 




y 2 


> = • 


^2 > 


0 


0 


M 3 








0 


^21 


^21 








0 


-C 
^21 


k 3 + k 2l _ 








l F 3, 



The trivial solution of the equation can be written as: 

AY" + BY' + CY = F 

It becomes, after multiplication by A" : 

Y" + A 1 BY' + A l CY = A l F 

We work in the case without default so Fl = F2 = F3 = 0, 

Therefore Y" + A * BY' + A l CY = A l F 
A" 1 is calculated by the shape: 



A-'=- 



1 



-coA' 



And A 1 = 



detA 

With: 

det(A) = M l M 1 M 3 

0 

0 
1 

M 3 

By replacing matrices by actual experimental values of the masses (Mi, M 2 , M 3 ), damping (ki, k 3 , k 2 i, k 2 3) and 

stiffness (C 2 i, C 23 ). 

ki, k 2 calculated from (5) 

C21, C23: calculated from (6) 

Mi = 25*10" 3 kg, M 2 = 2*10" 3 kg, M 3 = 23*10" 3 Kg 

k, = 10 7 N/m, k 3 = 10 7 N/m, k 2] = 0, 9925 N/m, k 23 = 0, 9925 N/m 

C 21 = 9, 9072, C 23 = 9, 9072 



1 






0 


m[ 






1 


0 




M 2 


0 


0 
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The mass matrix becomes: 

r 0 0 0 



B = 



With: 



0 9,9 -9,9 
0 -9,9 9,9 



A l B = 



40 0 



A~ 1 B = 



0 

0 500 0 

v 0 0 43,5^ 

And: 

r 0 0 0 

0 4950 -4950 

v 0 -413,25 413,25 

The damping matrix: 
^10 7 0 
0 0,99 
0 -0,99 



0 


0 


0 


0 


9,9 


-9,9 


0 


-9,9 


9,9 



c 



0 

-0,99 



10 s 



J 



With: 



A l C 



And: 



^40 0 
0 500 



0 
0 



0 0 43,5 



10 7 
0 
0 



0 
0,99 
-0,99 



0 

-0,99 
10 8 



Uoao 7 0 0 

0 495 -495 
0 -43 435. 10 7 

V 

The matrix L which explicitly interpolating numerical application: 



A l C = 



L = 



( 0 


0 


0 


1 


0 


0 ^ 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


1 


40. 10 7 


0 


0 


0 


0 


0 


0 


495 


-495 


0 


4950 


-4950 


v 0 


-43 


435. 10 7 


0 


-413,25 


413,25 y 



In seeking solutions of X: 



A-I 6 X = 



The polynomial to be solved is given by: 



' -X 0 


0 


1 0 


0 


0 -X 
0 0 


0 

-X 


0 1 
0 0 


0 
1 


40. 10 7 0 


0 


-X 0 


0 


0 495 


-495 


0 4950 -X 


-4950 


0 -43 


435.10 7 


0 -413,25 


413,25- 
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P(X) = X 6 -5,363.X 5 -4,75.10 9 .X 4 + 2368.10 I3 .X 3 +1,74.10 I8 .X 2 -8,613.10 21 .X -8,61 
The own values are: 

V, =2.10 4 
V 2 =-2.10 4 
V 3 =-6,5763.10 4 
V 4 =6,6178. 10 4 
V 5 =4,948. 10 3 

v 6 = -o,i 

The vectors associated with each own values are: 

^ = 5.10" 5 ^ + ^ 

v 1 =5.w i 7 l -7 i 

^ = 1,0618.10 ^ + 1,517.10 -6,98.10 2 ^ -9,97.10 ^ 
V 4 = -1, 22. 10 - 6 e~l - 1, 5. 10 " 5 ^ + 8, 06. 1 0" 2 - 9, 97. 1 0 ~ l 7 6 
V 5 =2, 02. 10" 4 7 2 + 9, 55.10 8 7 3 + 7 5 + 4, 72. 10" 4 ^ 
t£ = -9,95.10 5 ^ + 3, 83. 1 0 10 - 9, 95. 10 2 -3,83. 10""^ 

The explicit matrix vectors calculated: 

'5.10 5 5.10 5 0 0 0 0 N 

0 0 1,0618.10 6 1, 22.10 6 2,02.10" 4 9,95.10"' 

0 0 1,517.10 5 -1,5.10 5 9,55.10 8 3,83.10 10 

1 -1 0 0 0 0 

0 0 -6,98.10 2 8,06.10 2 1 -9,95.10 2 
v 0 0 -9,97.10"' -9,97.10"' 4J2.10" 4 -3,83.10 n y 

With: 

P: matrix that contains the own vector 
Furthermore P" 1 results in the following matrix: 



10 4 


0 


0 


85.10 1 


0 


0 ) 


10 4 


0 


0 


-5.10 1 


0 


0 


0 


-3,05.10 4 


3,31. 10 4 


0 


-2,93.10 3 


-0,5 


0 


-3,53.10 4 


-3,31. 10 4 


0 


3,4.10 3 


-0,5 


0 


9,994.10 2 


4,98.10 3 


0 


1 


5,83.10 3 


v0 


1 


-1 


0 


-2,02.10" 4 


-3,47.10" 8 y 



And: 
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rAP-- 



2.10 4 0 0 0 

0 -2.10 4 0 0 

0 0 -6,5763.10 4 0 

0 0 0 6,6178.10 4 



4948.10 3 



-0,lj 



And we consider : 



V 



y 2 
y'i 
-y\ 
-y's 
-y 6 



This vector can be written as follows: 

f 



y'i 


0 


0 


0 


0 


o N 


0 


y'i 


0 


0 


0 


0 


0 


0 


^3 


0 


0 


0 


0 


0 


0 


-yl 


0 


0 


0 


0 


0 


0 


-y'i 


0 


0 


0 


0 


0 


0 


-y'u 



= e v '+K 



With the vector V can be written as: 



v = 



\ 



0 

A l C A l B 



The previous system is reflected in the resolution of the matrix: 



0 j, 0 0 0 0 

0 0 j, 0 0 0 

0 0 0 -jj 0 0 

0 0 0 0 -j, 0 

0 0 0 0 0 -y h 

y[ = e 2AQit +1 
Incorporating includes: 

2.10 4 J 



„-6,5763.10 4 l 



tl 



tl 



Hence the withdraws: 



yi = 



2.10 4 



■ + t + A 



With the initial conditions imposed on the constant A to the value: 

1 

y 1= 0 A = — 



At t = 0 then 

Hence the solution 



2.10 4 
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e 2A0t 1 

y, = T + t r (10) 

^ 2.10 4 2.10 4 

Similarly for yl we find the general expression of y 2 : 



-2.Vft j 



g -2.10S 



y 7 = + 1 + B 

2 2.10 4 

With the initial conditions imposed there is the value of the constant B: 
At t = 0 y2=0 B 1 



2.10 4 



Where: 



e- 2AQit I 



y? = r +t + r (11) 

2.10 4 2.10 4 

Integrating is the solution of the vector y 3 : 
^=^ 57631 ° 4 <+1 

e -6,5763.10 4 r 

= 7 +t + C 

6,5763.10 4 

With the initial conditions imposed there is the value of the constant C: 
At t = 0 V, = 0 C = 

6,5763.10 4 



-6,5763. 10 4 f 



1 



y " 6,5763.10 4 +? 6,5763.10 4 

V. CONCLUSION 

Vibration analysis, for conditional preventive maintenance, proves a wonderful tool for several decades 
for the industry. It serves three levels of analysis: monitoring, diagnosis and damage of equipment state. 
This work helped to develop an analytical model that simulates the behavior of a vibratory ball bearing with 
three degrees of freedom.The problem addressed in this paper is the development of an analytical simulator 
vibration bearing. Theoretical research is detailed enough to provide a simple model that lends itself easily to 
programming, but powerful enough to incorporate the effect of a maximum of influential parameters on the 
vibration of bearings. 
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